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Abstract 

Vbfnlo is a flexible parton level Monte Carlo program for the simulation of 
vector boson fusion, double and triple vector boson production in hadronic collisions 
at next-to-leading order (NLO) in the strong coupling constant, as well as Higgs 
boson plus two jet production via gluon fusion at the one-loop level. In the new 
release - VERSION 2.5.0 - several new processes have been added at NLO QCD: 
vector boson fusion production of a Higgs boson plus a photon, vector boson fusion 
production of a photon, W^ and WZ production plus a hadronic jet and the triboson 
production processes WWj, ZZj, WZj, W^-y, Zjj and 777. The code has been 
extended to run in the Minimal Supersymmetric Standard Model (MSSM), and 
electroweak corrections to Higgs boson production via weak boson fusion have been 
included. Anomalous gauge boson couplings can be used in new processes and 
the Three-Site Higgsless model has been implemented for several processes. The 
simulation of Higgs boson production via gluon fusion has been improved. 
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1 Introduction 

The physics potential of the TeVatron and even more of the Lhc depends heavily on 
our ability to provide accurate cross section predictions for both signal and background 
processes. The latter are often generated by QCD interactions which give rise to weak 
bosons in the final state. A precise description of such hard QCD production processes is 
needed, as well as a method for simulating the measurable hadronic final states. Reaching 
these goals requires next-to-leading order (NLO) QCD calculations presented in the form 
of parton level Monte Carlo (MC) generators, which are an efficient solution when it 
comes to final states characterized by a high number of jets and/or identified particles. 
When kinematic cuts are imposed, as is mandatory for processes involving QCD radiation, 
analytic phase-space integration becomes impractical and implementation of results in the 
form of Monte Carlo programs becomes the method of choice. 

Vbfnlo is a fiexible MC program for vector boson fusion (VBF), double and triple 
vector boson production processes at NLO QCD accuracy. Furthermore, the electroweak 
corrections to Higgs boson production via VBF (which are of the same order of magni- 
tude as the QCD corrections in the experimentally accessible regions of phase-space) have 
been included. Since real emission processes are part of the NLO cross sections, Vbfnlo 
provides the means to calculate cross sections for the corresponding process with one ad- 
ditional jet at leading order (LO) in the strong coupling. In addition, the simulation of 
CV-even and CP-odd Higgs boson production in gluon fusion, associated with two addi- 
tional jets, is implemented at LO QCD. The full top- and bottom-quark mass dependence 
of the one- loop contributions in the Standard Model (SM), in the Minimal Supersym- 
metric Standard Model (MSSM) and in a generic two-Higgs-doublet model is included. 
Vbfnlo can be run in the MSSM (with real or complex parameters), and anomalous cou- 
plings of the Higgs boson and gauge bosons have been implemented for certain processes. 
Additionally, two Higgsless extra dimension models are included - the Warped Higgsless 
scenario and a Three-Site Higgsless Model - for selected processes. These models can be 
used to simulate the production of technicolor-type vector resonances in VBF and triple 
weak boson production. 

Arbitrary cuts can be specified as well as various scale choices. Any currently available 
parton distribution function (PDF) set can be used through the Lhapdf library. In 
addition, CTEQ6L1 for LO and CTIO for NLO calculations, as well as MRST2004qed 
and MSTW2008, are hard-wired into the code. For most processes implemented at leading 
order the program is capable of generating event files in the Les Houches Accord (LHA) 
and the HepMC format. When working in the MSSM, the SUSY parameters can be input 
via a standard SLHA file. 

This manual supersedes the previous version [1] released with Vbfnlo VERSION 2.0. 
The following additional processes have been included to NLO QCD accuracy: 

• Vector boson fusion production of a Higgs boson in association with a photon and 
two jets, 

• VBF production of a photon plus two jets, 

• diboson (WZ and W'y) plus a hadronic jet production, 

• WW'-f, ZZ'-f, WZ'-f triboson production, 

• W^'-ff, Z'-f'-f triboson production. 



• 777 production. 

Additionally, various new features have been introduced: 

• Anomalous quartic gauge boson couplings for triboson production, 

• Extension to the MSSM, 

• Electroweak corrections to VBF Higgs boson production in the SM and MSSM, 

• Triple vector boson production in Higgsless Kaluza-Klein models. 



The Vbfnlo webpage - http://www-itp.particle.uni-karlsruhe.de/vbfnlo/ 



contains, in addition to the latest version of the code, extra information such as the explicit 
implementation of the electroweak parameters and couplings, as well as optimised grid 
files for all processes for a set of standard cuts. To enable a simple installation test 
Vbfnlo is shipped with a complete set of example results, together with input files, in 
the regress directory. 



2 Installing Vbfnlo 

The source code of the current version of Vbfnlo can be downloaded from the Vbfnlo 
webpage 



http : //www-itp . particle . uni-karlsruhe . de/vbf nlo/ 



and includes a GNU conforming build system for portability and an easy build and 
installation procedure. 

2.1 Prerequisites 

The basic installation requires GNU make a FORTRAN770and a C++ compiler. Vbfnlo 
offers the possibility of using the Lhapdfo P] library for parton distribution functions. In 
order to include the electroweak corrections, the program LOOPTOOL^J |21ll] is required. 
Additionally, FeynHigg^j jlHH] can be linked to the code in order to calculate the Higgs 
boson sector of the MSSM, although a SLHA file can be used as an alternative. If the 
simulation of Kaluza-Klein resonances is enabled, an installation of the GNU Scientific Li- 
brary (GSL J3 is required. Vbfnlo can also be linked to ROOTO and HepMC0 to produce 
histograms and event files in those formats. 

2.2 Compilation and installation 

After unpacking the source archive and entering the source directory, the configure script 
can be invoked with several options, a complete list of which are available via . /configure 
--help. Among these, the most important ones are: 

• - -pref ix= [path] 

Install Vbfnlo in the location given by [path] . 

• --enable-processes=[list] 

By default, the code for all available processes is compiled. Optionally, [list] gives 
a comma-separated list of selected processes to be compiled. Possible process names 
are: 
vbf Vector boson fusion processes 

diboson Double gauge boson production 

triboson Triple gauge boson production 

diboson jet Double gauge boson production with a hadronic jet 
hj j j Higgs boson plus three jet production via vector boson fusion 

ggf Higgs boson plus two jet production via gluon fusion 

• --disable-NLO 

Disable the next-to-leading order QCD corrections. 



^The following compilers have been tested; g77 and gfortran. The ifort compiler has also been 
tested with a basic installation (i.e. without linking to any external programs). 
" http : //pro j ects . hepf orge . org/lhapdf/l 
http : //www . f eynarts . de/looptools/i 



http : //www . f eynhiggs . de/ 



http : //www . gnu . org/sof tware/gsl/ 



http : //root . cern . ch/ 



http: //Icgapp. cern. ch/project/simu/HepMC/ 



• --enable-kk 

Enable simulation of Kaluza-Klein resonances. Disabled by default, the Kaluza- 
Klein option requires the installation of the GNU Scientific Library, which needs to 
be specified via --with-gsl. 

• --with-gsl=[path] 

Enable the use of the GNU Scientific Library, [path] specifies the location of the 
GSL installation. 

• --with-LHAPDF=[path] 

Enable the use of Lhapdf in addition to the built-in PDF sets. Disabled by default, 
[path] specifies the location of the Lhapdf installation. 

• --with-LOOPTOOLS=[path] 

Enable the use of LOOPTOOLS in order to calculate the electroweak corrections. If 
this option is not specified, the electroweak corrections cannot be included. Disabled 
by default, [path] specifies the location of the LoOPToOLS installation. 

• --with-FEYNHIGGS=[path] 

Enable the use of FeynHigg^j to calculate the MSSM Higgs sector parameters. 
Disabled by default, [path] specifies the location of the FeynHiggs installation. 

• --with-root=[path] 

Enable the use of ROOT, [path] specifies the location of the ROOT installation. 

• --with-hepmc=[path] 

Enable the production of HepMC format event files, [path] specifies the location 
of the HepMC installation. 

Note that, by default, both LoOPToOLS and FeynHiggs are installed as static li- 
braries. If this is the case, configure must be run with the option --enable-shared=no. 
Also note that, in order to link to an external program such as LoOPToOLS, the exter- 
nal program needs to have been compiled using the same compiler (e.g. gf ortran) as 
Vbfnlo. Once configure has finished successfully, the make and make install com- 
mands will compile and install Vbfnlo, respectively. 

2.3 Source and installation directory layout 

The Vbfnlo source tree contains the following subdirectories: 

• amplitudes/: Routines to calculate matrix elements for the processes provided. 

• doc/ : The source of this manual. 

• helas/ : Helas j9] subroutines used to calculate helicity amplitudes. 

• loops/: One- loop tensor integrals up to five-point functions |1U| . 

• PDFsets/: Built-in parton distributions (CTEQ6L1 [II] for LO and CTIO [I2] for 
NLO calculations, as weU as MRST2004qed [I3] and MSTW2008 jUj). 



^Note that different versions of FeynHiggs give slightly different results as more corrections are added 
to the calculations. Vbfnlo has been tested with FeynHiggs versions 2.6.x, 2.7.x, 2.8.0 and 2.8.3. 



• phasespace/: Specialized phase-space generators for the processes provided. 

• regress/: Folder containing example results, together with input files, for all pro- 
cesses. 

• src/ : Source code of the main programs and input files. 

• utilities/: Routines for administrative tasks, cuts, scale choices and interfaces. 

The source does not need to be modified to change the simulation parameters. Vbfnlo 
offers several kinematic cuts and scale choices. This is illustrated in Sec. HI In ad- 
dition, it provides a few basic histograms. Cuts, histograms and scale choices not al- 
ready provided may be added in the utilities/cuts. F, utilities/histograms. F and 
utilities/scales. F files. 

The installation is performed in a standard UNIX-layout, i.e. the directory specified 
with the --prefix option of the configure script contains the following directories: 

• bin/: vbfnlo and ggflo executables. 

• include/VBFNLO/ : Vbfnlo header files. 

• lib/VBFNLO/: Vbfnlo modules as dynamically loadable libraries. These can also 
be used independently from one of the main programs. 

• share/VBFNLO/ : Input files and internal PDF tables. 

2.4 Running the program 

Both the vbfnlo and ggflo executables contained in the bin directory of the installation 
path look for input files in their current working directory. An alternative path to input 
files may be specified explicitly by passing the --input= [path] argument to the programs, 
with path denoting the full path where input files are located. I.e. in order to run Vbfnlo 
from the installation (prefix) directory, the command is 

./bin/vbfnlo --input=[path] 

The input files contained in the share/VBFNLO directory are meant to represent default 
settings and should not be changed. We therefore recommend that the user copies the 
input files to a separate directory. Here, special settings may be chosen in the input files 
and the program can be run in that directory without specifying further options. 

Vbfnlo outputs a running 'log' to the terminal, containing information about the 
settings used. In addition, the file xsection.out is produced, which contains only the 
LO and NLO cross sections, with the associated errors. Histograms and event files, in 
various forms, can be output as described later. 

2.5 MacOSX 

Unfortunately, owing to the library set-up, at the moment Vbfnlo does not run on 
MacOSX. It has, however, been tested successfully on a virtual box on a MacOSX. 



2.6 Bug reports 

Please report any problems to 

vbf nloOparticle . uni-karlsruhe . de 

with a short report including the configure options used to build Vbfnlo, as well as the 
versions of compilers and external libraries used. 

2.7 License 

Vbfnlo is distributed under the GNU General Public License (GPL) version 2. This 
ensures that the source code will be available to users, grants them the freedom to use 
and modify the program and sets out the conditions under which it can be redistributed. 
However, it was developed as part of an academic research project and is the result of 
many years of work by the authors, which raises various issues that are not covered by the 
legal framework of the GPL. It is therefore distributed together with a set of guidelineqj, 
which were originally formulated and agreed on by the MCnet collaboration for event 
generator software. 



='These guidelines are contained in the GUIDELINES file distributed with the release. 



3 Processes 



In the following sections, we describe all production processes and decay modes imple- 
mented in Vbfnlo, together with references to more detailed discussions of the underlying 
calculations. 

In the phase-space regions that are accessible at hadron colliders, VBF reactions are 
dominated by t-channel electroweak gauge boson exchange. In Vbfnlo s-channel ex- 
change contributions and kinematically suppressed fermion interference contributions [T5] - 
[T7] are therefore disregarded. Numerically small contributions from Pauli-interference ef- 
fects for identical charged leptons are neglected. For final-state identified photons we 
employ the isolation criterion of Ref. |18] . This ensures that divergences from collinear 
emission of a photon from a final-state massless quark or anti-quark are avoided, while 
the cancellation of the infrared divergences between the virtual and the real part is kept 
intact. 



3.1 VBF Higgs boson production in association with two jets 

Hjj production via VBF mainly proceeds via electroweak quark-quark scattering pro- 
cesses such as qq' — > qq'H and crossing related reactions. In Vbfnlo, tree level cross 
sections, NLO QCD corrections and 1-loop electroweak corrections (in the Standard Model 
and the MSSM) to the t-channel production process are provided. The subsequent de- 
cay of the Higgs boson is simulated in the narrow width approximation (NWA). For the 
H -^ W^W~ and the H — )■ ZZ modes, full off-shell effects and spin correlations of the 
decay leptons are included. The available production process and decay modes are listed 
with the corresponding process IDs in Table [1] Anomalous couplings between a Higgs 
boson and a pair of vector bosons are implemented in the code. Details of the calculations 
can be found in Refs. ll9f[2T]. 



ProcId 


Process 


BSM 


100 


PP -^ Hjj 




>! 






101 


PP -^ Hjj 


^11 jj 








102 


PP -^ Hjj 


-^ f^^fJ-' 33 








103 
104 


PP -^ Hjj 
PP -^ Hjj 


-^ T+r~ jj 
-^ bbjj 




> anomalous HVV couplings, 


MSSM 


105 


PP -^ Hjj 


-^W+W~ jj^ it ^h^2^i2 3 3 








106 


PP -^ Hjj 


-^zzjj^eii^e+qjj 








107 


PP -^ Hjj 


-^ ZZjj -^ i^i^i^i^Pe^ jj 


- 







(-) 



Table 1: Process IDs for pp —^ Hjj production via weak boson fusion at NLO (QCD 
and electroweak) accuracy in the SM and MSSM. Anomalous Higgs boson couplings are 
implemented for all decay modes. 



ProcId 


Process 


110 


(-) 


^Hjjj 




111 


(-) 
PP 


-^Hjjj 


-^ 11333 


112 


(-) 
PP 


^Hjjj 


-^ fJ-^fJ-'JJJ 


113 


(-) 
PP 


-^Hjjj 


-> T+r~ jjj 


114 


(-) 
PP 


-^Hjjj 


-^ bbjjj 


115 


(-) 
PP 


-^Hjjj 


-^ W+W-jjj -^ It^ej2^i2 333 


116 


{-) 
PP 


^Hjjj 


^zzjjj^i+i^itqjjj 


117 


(-) 
PP 


-^Hjjj 


-^ zzjjj -^ et^i^h'^hnj 



(-) 



Table 2: Process IDs for pp — t- Hjjj production via VBF at NLO QCD accuracy. 



ProcId 


Process 


2100 


PP -^ H-fjj 




2101 


PP -^ H-fjj 


-^ 11133 


2102 


PP ^ H-fjj 


^ ^-^fi^ljj 


2103 


PP -^ H-fjj 


-^T+T-'jjj 


2104 


PP -^ H-ijj 


-^ bb-fjj 


2105 


PP -^ H-fjj 


-^ w+w-jjj -^ iti^ij2^i2in 


2106 


PP ^ H-fjj 


-^ ZZ^jj ^ it 1^4^133 


2107 


PP ^ Hjjj 


-^ zZ'jjj -^ elqi^i^Pi^-fjj 



(-) 



Table 3: Process IDs for pp — )■ H'jjj production via VBF at NLO QCD accuracy. 



3.2 VBF Higgs boson production in association with three jets 

Adding an extra parton to the Higgs production processes of Sec. 13.11 gives rise to Hjjj 
final states. The corresponding cross sections are implemented at NLO QCD accuracy 
(with certain approximations) in Vbfnlo. A list of all available modes and corresponding 
process IDs is given in Table |2j Details of the calculation can be found in Ref. |22) . 



3.3 VBF Higgs boson production with a photon and two jets 



The emission of an additional photon in VBF Higgs boson production (Sec. 13. ip results 
in H'-fjj final states. These are implemented at NLO QCD accuracy in Vbfnlo, with 
process IDs as given in Table [3l Details of the calculation can be found in Ref. |23] . 



ProcId 


Process 


120 


{-) 
PP 


^zjj^e+£-jj 


121 


(-) 
PP 


-^ Zjj -^ vivaj 


130 


{-) 
PP 


-^w+jj^£+uaj 


140 


(-) 
PP 


-^W" jj -^(.-Viij 


150 


{-) 
PP 


-^in 



Table 4: Process IDs for vector boson 
QCD accuracy. 



2 jet production via weak boson fusion at NLO 



3.4 VBF production of a vector boson and two jets 

Vector boson fusion processes can also produce final states with two leptons plus two jets, 
which are generically referred to as "VBF Zjj and W^jj production". These reactions 
and the one with a photon plus two jets in the final state are implemented to NLO QCD 
accuracy in Vbfnlo, see Table |H Details of the calculations can be found in Refs. [211125] ■ 



3.5 VBF production of two vector bosons and two jets 

The production of four leptons plus two jets in the final state at order 0{a^) is dominated 
by VBF contributions. In Vbfnlo, all resonant and non- resonant t-channel exchange 
contributions giving rise to a specific leptonic final state are considered. For simplicity, 
we refer to these reactions as "VBF diboson production". Finite width effects of the 
weak bosons and spin correlations of the decay leptons are fully retained. The available 
processes and corresponding process IDs are listed in Table O Anomalous gauge boson 
couplings are implemented for W^W~ jj production, and this class of processes can also 
be run in the Higgsless models implemented in Vbfnlo. Details of the calculations can 
be found in Refs. [26H28]. 



3.6 Double and triple vector boson production 

The production of four- and six-lepton final states mainly proceeds via double and triple 
vector boson production with subsequent decays. Additionally, there are processes where 
one or more decaying boson is replaced by an on-shell photon, giving rise to lepton pro- 
duction in association with photon(s) and triple photon production. In Vbfnlo, the 
processes listed in Tableware implemented to NLO QCD accuracy, including full off-shell 
effects and spin correlations of the final state leptons and photons. For processes with 
three massive gauge bosons, the Higgs boson contributions are included in the Standard 
Modetlj. Anomalous vector boson couplings are implemented for W~^W~Z, ZZW^ and 



"'^^This should provide a good approximation in the MSSM in the decouphng regime - i.e. where one 
of the Higgs bosons is SM-hke. 
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ProcId 


Process 


BSM 


200 


pp -^ W+W- jj -^ it^eJ2^i2JJ 


anomalous couplings, Kaluza-Klein 


models 


210 


Pp'^ZZjj^iil^it^jj 


' 






211 

220 


pp ^ zzjj -^ it^i^h^t2n 

p''^^W+Zjj^itiyi,iti2JJ 




> Kaluza-Klein models 




230 


p'-j^^w~Zjj^qu,,it^2 3J 


' 







Table 5: Process IDs for dihoson + 2 jet production via weak boson fusion at NLO QCD 
accuracy. 



ProcId 


Process 


BSM 


300 


(-) 
PP - 


-> W+W- -^ ifuej^'^ii 




400 


(-) 
PP - 


-> W+W-Z -^ efiye/^i?ej3+e3- 


>! 




410 
420 
430 


(-) 
PP - 

(-) 
PP - 

(-) 
PP - 


-> ZZW+ -^ efi^i^^^h^^is 

-> ZZW~ -^ £ti^ip2^3~i^k 




anomalous couplings, 
Kaluza-Klein models 


440 


(-) 
PP - 


-> W~W+W- ^ £^i?e,£tuej3-u,. 


, 




450 


(-) 
PP - 


-> ZZZ -^ £^£+£^£+£3-13+ 




460 


(-) 
PP - 


-> W-W+j -^ £^i?i^£+ue^j 




470 


(-) 
PP - 


-> ZZ-y -^ £^£f£^£^-f 




480 


(-) 
PP - 


^ W+Zj ^ £tiye,i2 4l 




490 


(-) 
PP - 


-> W-Z-f -^ ^^%£^47 




500 


(-) 
PP - 


-)■ VF+77 — )■ ^~'~z/f77 




510 


(-) 
PP - 


-7- VF~77 —7- ^~z?£77 




520 


(-) 
PP - 


-)■ Z77 — ;> £~^+77 




521 


(-) 
PP - 


-7- ^77 —7- i^ii'i'j'y 




530 


(-) 
PP - 


-)■ 777 





Table 6: Process IDs for the dihoson and triboson production processes at NLO QCD 
accuracy. 



W^W~^W production. These processes can also be run in the implemented Kaluza-Klein 
models. Details of the calculations can be found in Refs. |29H3' 
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ProcId 


Process 


BSM 


610 


pp — 7- M/ 7J - 


-^ ^ j^nj 


>, 




620 
630 


PP — ;> W^^J - 
PP —7- W~Zj 






anomalous couplings 


640 


PP -^ W^Zj 


^ ^|zV£,^2"^^i 


■ 





Table 7: Process IDs for W'~fj and WZj production at NLO QCD. 



ProcId 


Process 


BSM 


4100 


{-) „ . . 
PP ^ Hjj 




MSSM, general 2HDM 


4101 


PP ^ Hjj ^-f-fjj 


• 




4102 


PP ^ Hjj - 


-> /^"^/"~ jj 






4103 


PP -^ Hjj - 


-> -r+T^ jj 






4104 


PP ^ Hjj - 


-> bbjj 




. MSSM 


4105 


PP -^ Hjj - 


^W+W-jj^ttvi.qyiJj 






4106 


PP ^ Hjj - 


^zzjj^et^i4^2n 






4107 


PP -^ Hjj - 


-> ZZjj -^ ifl^iyi^fye^jj 


■ 





Table 8: Process IDs for LO Higgs boson plus 2 jet production via gluon fusion. 



3.7 Double vector boson production in association with a hadronic 
jet 

14^7 and WZ production in association with a hard hadronic jet are available in Vbfnlo 
at NLO QCD accuracy under the process IDs of Table [71 All off-shell and finite width 
effects in the fixed width scheme are included. The processes can be run with anomalous 
WW Z and WW'-f couplings. Details can be found in Refs. 



3.8 Higgs boson production in gluon fusion with two jets 

CP-even and CP-odd Higgs boson production in gluon fusion, associated with two addi- 
tional jets, is a process that first appears at the 1-loop level which, therefore, is counted 
as leading order in the strong coupling. This process is simulated with the full mass 
dependence of the top- and bottom-quark running in the loop in the Standard Model, 
in the (complex) MSSM and in a generic two- Higgs- doublet model. The relevant process 
IDs are given in Table |8l Details of the calculations can be found in Refs. 
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4 Input files and parameters 

Vbfnlo is steered through the following input files: 

• vbfnlo . dat : General parameters for a run. 

• cuts. dat: Values for kinematic cuts. 

• ggflo.dat: Additional parameters for the ggflo program. 

• susy.dat: Parameters describing the MSSM scenario. 

• anom_HVV . dat : Parameters for anomalous Higgs boson couplings. 

• anomV . dat : Parameters for anomalous gauge boson couplings. 

• kk_input . dat : Settings for the Warped Higgsless and Three-Site Higgsless Models. 

• kk_coupl_inp.dat: Numerical values if externally calculated Kaluza-Klein cou- 
plings and masses should be used. 

• histograms . dat : Histogram ranges. 

• random . dat : Seed for the random number generator. 

The following subsections give a detailed description of all available parameters. 

4.1 vbfnlo.dat 

vbfnlo.dat is the main input file for Vbfnlo. 

4.1.1 vbfnlo.dat — general parameters 

• PROCESS: Process ID as described in Sec. [31 Default is 100: Higgs boson production 
via VBF. 

• LOPROCESS_PLUS_JET: If set to true, the leading order process with one additional 
jet is generated, i.e. only the real radiation contribution is generated. This option 
is available for all but gluon fusion processes. Default is false. 

• LEPTONS: Choice of the final state leptons (decay products of W and Z bosons) 
according to the MC particle numbering scheme 
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J^e 
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f^~ 


14 


V,. 


15 


T^ 


16 


Ur 



If the selected configuration is not available, the appropriate first-generation leptons 
(electrons and electron neutrinos) are used as default values. 
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LO_ITERATIONS : Sets the number of iterations for the integration of LO cross sec- 
tions. Usually more than one iteration is used in order to adapt the integration grid 
and thus improve the efficiency of the MC integration algorithno. For an adapted 
grid file (see LO_GRID) this parameter can be set to 1. Default is 4. 

NLO.ITERATIONS: Analogous to LO_ITERATIONS, but for the real emission part of 
an NLO calculation. Since the corresponding phase-space is different from the LO 
configuration, a second independent MC integration has to be performed. Default 
is 4. 

LO_POINTS: Determines the number of phase-space points that are generated in 
each iteration. In the last iteration there are 2'^ points, where N = LO_POINTS. In 
each previous iteration, the number of points is half the value of the following one. 
Example: For 4 iterations (LO.ITERATIONS = 4) and LO.POINTS = 20, there are 
2i7 generated points in the first, 2^^ in the second , 2^^ in the third and 2^° ^ 10^ 
in the last iteratioiri. Default is A^ = 20. 



NLO.POINTS: Similar to LO_POINTS, but for the real emission part of an NLO cal- 
culation. Default is 20. 

LO_GRID: Sets the name of the grid files that are generated at the end of each 
iteration. Choosing name as the input parameter, in each iteration X a grid file 
name . out . X will be produced. If a grid file name is already present in the working 
directory, the program reads in this file when executed. Note that optimised grids 
for all processes (using standard cuts) are provided on the Vbfnlo webpage. 

NLO_GRID : Similar to LO_GRID, but for the real emission part of an NLO calculation. 

PHTN_GRID : Similar to NLO_GRID, but for the real photon emission part of an NLO 
electroweak calculation. 

NLO_SWITCH: Switch for the NLO part of a process, if available. If set to true, 
cross sections and histograms are calculated to NLO QCD accuracy. Default is set 
to false. 

EWCOR_SWITCH : Switch for the electroweak corrections (note that this is only avail- 
able for VBF Higgs boson production). If set to true, cross sections and histograms 
are calculated to NLO electroweak accuracy. This option can only be used if LOOP- 
TOOLS was enabled at compilation. Default is set to false. 

ECM : The center-of-mass energy, ^/s, of the collider, measured in GeV. Default is 
14000 GeV. 

BEAMl, BEAM2: Define the type of particle of each beam. Possible options are -|-1 
for proton beams and —1 for anti-proton beams. Default is proton-proton collisions, 

(+1, +1). 

ID_MUF : Choice of the factorization scale. See Table IHlfor a list of available options. 
Default is 0. 



^^For all NLO calculations the virtual contributions are calculated using the already optimized leading 
order grid. 

^^The virtual contributions are calculated for 2^ points only. 
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ID.MUF 


Factorization Scale 


Process class 



1 
2 
3 
4 


user defined constant scale set by MUF_USER 
momentum transfer of exchanged W/Z boson 

min(pT(j*)) 

invariant mass of VV system 

invariant mass of VVV system 


all 

vbf (except Hjjj) 

vbf 

diboson 

triboson 


5 
6 

7 


VPt(Ji) X Pt(J2) 

constant scale = Higgs boson mass 
Minimum transverse energy of bosons 


ggf 
all 

diboson, triboson 



Table 9: Factorization scale options. 



ID.MUR 


Renormalization Scale 


Process class 



1 
2 
3 
4 
5 
6 
7 


user defined constant scale set by MUR_USER 

momentum transfer of exchanged W/Z boson 

mm{p^{ji)) 

invariant mass of VV system 

invariant mass of VVV system 

at = asipAJi)) X as{pT{J2)) x aUniH) 
constant scale = Higgs boson mass 
Minimum transverse energy of bosons 


all 

vbf (except Hjjj) 

vbf 

diboson 

triboson 

ggf 
all 

diboson, triboson 



Table 10: Renormalization scale options. 



• ID_MUR: Choice of the renormalization scale. See Table [10] for a list of available 
options. Default is 0. 

• MUF_USER: If ID_MUF is set to 0, this parameter sets the user defined constant 
factorization scale measured in GeV. Default is 100 GeV. 

• MUR_USER: If 1D_MUR is set to 0, this parameter sets the user defined constant 
renormalization scale measured in GeV. Default is 100 GeV. 

• XIF: Factor by which the factorization scale is multiplied. May be used to analyze 
the scale dependence of differential cross sections. Default is 1. 

• XIR: Factor by which the renormalization scale is multiplied. May be used to 
analyze the scale dependence of differential cross sections. Default is 1. 

Note that alternative scale choices can be implemented in the file utilities/scales. F. 
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4.1.2 vbfnlo.dat — physics parameters 

• HMASS : Standard Model Higgs boson mass in GeV. Default value is 120 GeV. 

• HTYPE : Type of Higgs boson produced: 

- HTYPE = : SM Higgs boson, with mass HMASS 

- HTYPE = 1 : Light CP-even MSSM Higgs boson 

- HTYPE = 2 : Heavy CP-even MSSM Higgs boson 

- HTYPE = 3 : CP-odd MSSM Higgs boson (Note: this is not produced at LO 
unless we are working in the MSSM with complex parameters and include 
Higgs-propagator effects at LO.) 

Note that for HTYPE = 1-3, the input HMASS is not used. Default value is (SM 
Higgs boson). 

• MODEL: This flag determines whether we are working in the SM (1) or MSSM (2). 
Defauh is SM (1). Note that if HTYPE = 1-3 is chosen with MODEL = 1, the code 
will run in the Standard Model, but with a Higgs boson mass equal to that given 
by the specified MSSM parameters. 

• TOPMASS: Top-quark mass in GeV. Defauh value is 172.4 GeV. If a SLHA file is 
being used, the SLHA value will be taken rather than the input TOPMASS. 

• BOTTOMMASS : Bottom-quark pole mass in GeV, used in the calculation of the Higgs 
boson width and branching ratios. In the gluon fusion process, mb{MH) is used (cal- 
culated from the input pole mass). Default value is 4.855 GeV, which corresponds to 
m^^{mh) = 4.204 GeV. The explicit formula used is given on the Vbfnlo webpage. 
If a SLHA file is being used, the SLHA value will be taken rather than the input 
BOTTOMMASS. 

• CHARMMASS : Charm-quark pole mass in GeV used in the calculation of the Higgs 
boson width and branching ratios. Default value is 1.65 GeV, corresponding to 
mf^(mj = 1.273 GeV. If a SLHA file is being used, the SLHA value will be taken 
rather than the input CHARMMASS. 

• EWSCHEME : Sets the scheme for the calculation of electroweak parameters. A sum- 
mary of the six available options is given in Table [TT] The implementation of 
EWSCHEME has been simplified in the current version of VbfnlO: full details of all 
changes and their effects, together with the explicit formulae used, are available on 
the Vbfnlo webpage. Note that if EWSCHEME = 4 is chosen, all variables in Ta- 
ble [TT] are taken as inputs. As the parameters are not independent, this can lead 
to problems if the input values are not consistent. In this scheme, all photon cou- 
plings are set according to the input variable INVALFA and all other couplings are 
set according to FERMI _CONST. Note also that the choice of EWSCHEME can have a 
large effect on the relative size of the electroweak corrections, as the charge renor- 
malization depends on the way in which the electromagnetic coupling in the LO 
cross section is parametrized. Default value is 3. 

• FERMI_CONST: Fermi constant, used as input for the calculation of electroweak pa- 
rameters in EWSCHEME = 2 and 3. Defauh value is 1.16637 x 10"^ GeV"^. If a SLHA 
file is being used, the SLHA value will be taken rather than the input FERMI _CONST. 
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EWSCHEME 


Parameter 


Default Value 


Input/Calculated 




FERMI.CONST 


1.16637 X 10"5 GeV-2 


Input 




INVALFA 


128.944341122 


Input 


1 


SIN2W 


0.230990 


Calculated 




WMASS 


79.9654 GeV 


Calculated 




ZMASS 


91.1876 GeV 


Input 




FERMI.CONST 


1.16637 X 10-^ GeV"^ 


Input 




INVALFA 


132.340643024 


Calculated 


2 


SIN2W 


0.222646 


Input 




WMASS 


80.3980 GeV 


Calculated 




ZMASS 


91.1876 GeV 


Input 




FERMI.CONST 


1.16637 X 10-5 GeV-2 


Input 




INVALFA 


132.340705199 


Calculated 


3 


SIN2W 


0.222646 


Calculated 




WMASS 


80.3980 GeV 


Input 




ZMASS 


91.1876 GeV 


Input 




FERMI.CONST 


1.16637 X 10-5 GeV"^ 


Input 




INVALFA 


137.035999679 


Input 


4 


SIN2W 


0.222646 


Input 




WMASS 


80.3980 GeV 


Input 




ZMASS 


91.1876 GeV 


Input 




INVALFA (ZMASS) 


128.944341122 


Input 


5 


SIN2W 


0.222646 


Calculated 




WMASS 


80.3980 GeV 


Input 




ZMASS 


91.1876 GeV 


Input 




INVALFA (0) 


137.035999679 


Input 


6 


SIN2W 


0.222646 


Calculated 




WMASS 


80.398 GeV 


Input 




ZMASS 


91.1876 GeV 


Input 



Table 11: Electroweak input parameter schemes. 
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INVALFA: One over the fine structure constant, used as input for EWSCHEME = 1, 
4, 5 and 6. Witliin tlie otlier scliemes tliis parameter is calculated. The default 
value depends on the choice of EWSCHEME, as given in Table HH If EWSCHEME = 5 
is chosen, the value of a should be a{Mz\ whereas if EWSCHEME = 6 is chosen, the 
value of a should be a(0). In order to ensure backwards compatibility with previous 
versions of Vbfnlo, as an alternative ALFA, the fine structure constant, can be used 
as an input in vbfnlo.dat, which is read and used only if INVALFA is not present. 
If a SLHA file is being used, the SLHA value will be taken rather than the input 
INVALFA or ALFA. 



DEL. ALFA: Value of Aa, where 



.(M.,^^ (1) 



This is used as input for EWSCHEME = 6. Default value is 0.059047686. Note that 
this is only used for the electroweak corrections during the calculation of the charge 
renormalization constant. 

SIN2W: Sinus squared of the weak mixing angle. Used as input for EWSCHEME = 2 
and 4. Within the other schemes this parameter is calculated. Default input value 
is 0.222646. If a SLHA file is being used, the SLHA value will be taken rather than 
the input SIN2W. 

WMASS : W boson mass in GeV. This parameter is calculated in EWSCHEME = 1 and 
2. Default input value is 80.398 GeV. If a SLHA file is being used, the SLHA value 
will be taken rather than the input WMASS. 

ZMASS: Z boson mass in GeV. Default value is 91.1876 GeV. If a SLHA file is being 
used, the SLHA value will be taken rather than the input ZMASS. 

ANOM_CPL : If set to true, anomalous Higgs boson or gauge boson couplings are used 
if available for the selected process. These are available for the 

— if jj and W^W~jj production processes in VBF 

— triple vector boson production processes WWZ, W^W^W^ and ZZW^ 

— W'-fj and WZj production processes 

Default is set to false. 

KK_MOD: Option for the Warped Higgsless Model and Three-Site Higgsless Model. 
This is available for all VVjj production modes in VBF and for the triboson pro- 
cesses VV^W+W-, W+W'Z and W^ZZ. Default is set to false. Note that this 
needs to be enabled at compilation, using the configure option --enable-kk. 

EW_APPROX: Option controlling the electroweak corrections in Hjj production via 
VBF. 

— EW_APPROX = : No approximations involved. This option is not available 
when working in the MSSM. 

— EW.APPROX = 1 : Only top/bottom (and stop/sbottom in the MSSM) loops 
are calculated. 
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— EW_APPROX = 2 : All fermion (and sfermion in the MSSM) loops are calculated. 

— EW_APPROX = 3 : MSSM option - all SM-type (i.e. fermions, gauge and Higgs 
bosons) and sfermion loops are calculated. 

- EW.APPROX = 4 : MSSM option - all MSSM corrections to the Higgs boson 
vertex are calculated, together with all SM-type and sfermion corrections else- 
where. 

- EW.APPROX = 5 : MSSM option - all MSSM corrections to the Higgs boson 
vertex, the quark vertex and the vector boson self energies are calculated, to- 
gether with all SM-type boxes and pentagons (i.e. only chargino and neutralino 
box and pentagon diagrams are neglected). 

When working in the SM, EW.APPROX options 3-5 are equivalent to EW_APPROX = 0. 
Default is 5: full corrections in the SM, and the most complete available corrections 
in the MSSM. 

4.1.3 vbfnlo.dat — parameters for event output 

Vbfnlo generates parton level events according to the most recent Les Houches Ac- 
cord (LHA) format [IHI and in the HepMC format jU] for processes at leading order. 

• LHA_SWITCH: Switch on or off output of LHA event files. Default is set to false. 
Note that LHA event file output is not yet available for diboson plus jet processes. 

• LHA_FILE: Name of output LHA event file. Default is "event. Ihe". 

• HEPMC.SWITCH: Switch on or off output of HepMC event files. Default is set to 
false. Note that HepMC event file output is not yet available for diboson plus jet 
processes. 

• HEPMC.FILE: Name of output HepMC event file. Default is "event . hepmc". 

• UNWEIGHTING_SWITCH : Option for event weights. If set to true, events are un- 
weighted (event weight = +1). If set to false, events are weighted. Default is set 
to false. 

• PRENEVUNW: The number of events used in the last iteration in order to calcu- 
late/estimate the premaximal weight that is needed in the first step of the un- 
weighting procedure. Default is 1000. After all events are unweighted, the maximal 
weight is again calculated and a reweighting procedure is applied. 

• TAUMASS : Option to include the mass of r leptons in the event files. This option is 
only available for the vbf processes. Default is set to false. 

4.1.4 vbfnlo.dat — PDF parameters 

Vbfnlo can use either built-in parton distribution functions (PDFs) or the Lhapdf 
library. 

• PDF.SWITCH: Option to choose which PDFs are used. If set to 0, built-in PDFs 
CTEQ6L1 [n] for LO and CTIO pS] for NLO calculations are used. If set to 2, 
built-in PDFs MRST2004qed [13] are used at LO and NLO (if this option is chosen, 
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photon-induced processes can be included when calculating the electroweak correc- 
tions to Hjj production via VBF) and if set to 3 the built-in PDFs MSTW2008 [H] 
are used. If set to 1, an interface to Lhapdf is provided. Default is 0. 

The following options are used if Lhapdf has been selected (i.e. PDF_SWITCH = 1). 

• LO_PDFNAME: Name of the LO PDF set, see PDFsets. index or Ref. j2]. Default is 
"cteqSll.LHpdf" (CTEQ6L1). 

• NLO.FDFNAME: Name of the NLO PDF set, see PDFsets. index or Ref. [2]. Defauh 
is "CTlO.LHgrid" (CTIO). 

• LO_PDFMEMBER: Member PDF of the LO PDF set. Default is 0. 

• NLO.PDFMEMBER: Member PDF of the NLO PDF set. Default is 0. 

For compatibility with earlier versions of Vbfnlo, the following two variables are also 
supported for setting the PDF sets when using Lhapdf. If both are present, LO_PDFNAME 
and NLO.PDFNAME take precedence over LO.PDFSET and NLO_PDFSET, respectively. 

• LO_PDFSET: LHAGLUE number for the LO PDF set, see PDFsets . index or Ref. [2]. 

• NLO.PDFSET: LHAGLUE number for the NLO PDF set, see PDFsets. index or 
Ref. 12]. 

4.1.5 vbfnlo.dat — parameters for histograms 

Vbfnlo provides output for histograms in the following formats: TOPDRAWEH^i. ROOT 
and Gnuplottj, as well as raw data tables. 



ROOT: Enable output of histograms in ROOT format. Default is set to false. Addi- 
tionally, custom Root histograms can be defined by the user in utilities/rootuserhisto . cpp, 
which are filled with the weighted events from Monte Carlo integration. This option 
needs to be enabled when building Vbfnlo using the option --with-root. 

TOP: Enable output of histograms in TOPDRAWER format. Default is set to false. 

GNU: Enable output of histograms in Gnuplot format. Default is set to true. 

DATA: Enable output of raw data in a directory hierarchy. Default is set to true. 

REPLACE: Switch to overwrite existing histogram output files. Default is set to true. 

ROOTFILE: Name of the ROOT output file. Default is histograms. 

TOPFILE: Name of the TOPDRAWER output file. Default is histograms. 

GNUFILE: Name of the Gnuplot output file. Default is histograms. 

DATAFILE: Name of the data output directory. Default is histograms. 



http : //www ■ pa . msu . edu/ref erence/topdrawer-docs/| 



http : //www . gnuplot . info/ 
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4.2 cuts.dat — parameters for kinematic cuts 

The following general set of cuts has been implemented in Vbfnlo. Alternative cuts can 
be added in the file utilities/cuts. F. 

Jet-specific cuts: 



• RJJ_MIN: Minimum separation of two identified jets, ARjj = w A|/?- + A(/)^-, used 

by the generalised k± jet finding algorithm [IS] that combines all partons. Default 
is 0.8. 

• Y_P_MAX: Maximum allowed pseudorapidity for observation of final state partons 
(detector edge). Default is 5.0. 

• PGENKT JET : Exponent of the generalised k± algorithm. This yields the k± algorithm 
when setting the variable to 1, the Cambridge/ Aachen algorithm [12] for and the 
anti-/cx algorithm [5U] when setting it to —1. Default value of the floating-point 
number is 1.0. 

• PT_JET_MIN: List of minimum transverse momenta for identified jets in descending 
order. The later values can be omitted if they are the same as previous values. 
Default is 20 GeV. 

• Y_JET_MAX: Maximum allowed rapidity for identified jets. Default is 4.5. 

Lepton specific cuts: 

• Y_L_MAX: Maximum pseudorapidity for charged leptons. Default is 2.5. 

• PT_L_MIN : Minimum transverse momentum for charged leptons. Default is 10 GeV. 

• MLL_MIN : Minimum invariant mass for any combination of oppositely charged lep- 
tons. Default is 15 GeV. 

• MLL_MAX : Maximum invariant mass for any combination of oppositely charged lep- 
tons. Default is 10^° GeV. 

• RLL_MIN: Minimum separation of charged lepton pairs, Ai?«. Default is 0.4. 

• RLL_MAX: Maximum separation of charged lepton pairs, ARu. Default is 50. 
Photon specific cuts: 

• Y_G_MAX: Maximum pseudorapidity for photons. Default is 1.5. 

• PT_G_MIN : Minimum transverse momentum for photons. Default is 20 GeV. 

• RGG_MIN: Minimum separation of photon pairs, AR^^. Default is 0.6. 

• RGG_MAX: Maximum separation of photon pairs, AR^^. Default is 50. 

• PHISOLCUT: Photon isolation cut 6o as defined in Ref. [H]. Defauh is 0.7. 
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• EFISOLCUT: Efficiency e of photon isolation cut. Default is 1. The photon isolation 
cut is then described by: 

V EtAS - Ri-y) < e PT-, ^—^ for all 5 < So (2) 

■^-^ 1 — cos do 

where i is a parton with transverse energy Et^ and a separation Ri^ with a photon 
of transverse momentum prp^. 

Additional cuts: 

• RJL_MIN: Minimum separation of an identified jet and a charged lepton, ARj£. 
Default is 0.6. 

• RJG_MIN: Minimum separation of an identified jet and a photon, ARj^. Default is 
0.6. 

• RLG_MIN: Minimum separation of a charged lepton and a photon, AR^^. Default is 
0.6. 

Jet veto: 

• JVETO: If set to true, a jet veto is applied. For vbf and ggf processes it is applied 
to central jets beyond the two tagging jets. For all other processes the jet veto is 
applied to additional jets beyond the leading-order number, ordered by decreasing 
transverse momentum. Default is false. 

• YMAX_VETO: Maximum pseudorapidity of the additional jet. Default is 4.5. 

• PTMIN_VETO: Minimum transverse momentum of the additional jet. Default is 
50 GeV. 

• DELY_ JVETO: Minimum rapidity separation of a central jet from the two tagging 
jets for vbf and ggf processes. Default is 0. 

VBF specific cutO 

• ETAJJ_MIN: Minimum required pseudorapidity gap, Arjjj, between two tagging jets 
(the two leading jets in a p^ ordering). Default is 0. 

• YSIGN: If set to true, the two tagging jets are required to be found in opposite 
detector hemispheres. Default is false. 

• LRAPIDGAP : If set to true all charged leptons are required to lie between the two 
tagging jets in rapidity. Default is false. 

• DELY_JL: Minimum rapidity distance of the charged leptons from the tagging jets, 
if LRAPIDGAP is set to true. Default is 0. 

• GRAPIDGAP : If set to true all photons are required to lie between the two tagging 
jets in rapidity. Default is false. 

• DELY_JG: Minimum rapidity distance of photons from tagging jets, if GRAPIDGAP is 
set to true. Default is 0. 

• MDIJ_MIN: Minimum dijet invariant mass of two tagging jets. Default is 600 GeV. 

• MDIJ_MAX: Maximum dijet invariant mass of two tagging jets. Default is 10^° GeV. 



^'"'These apply only to vbf and ggf processes. 
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4.3 ggflo.dat — general parameters for gluon fusion processes 

In Vbfnlo the double real-emission corrections to gg — )■ 0, which lead to + 2 jet events 
at order a^, are included |11]. Here, can be set to be a Standard Model Higgs boson 
or any of the three neutral MSSM Higgs bosons (including mixing between h, H and A 
in the real or complex MSSM) by using the variables HTYPE and MODEL in vbfnlo.dat, 
or a mixture of scalar and pseudoscalar Higgs bosons as in a generic two-Higgs-doublet 
model (2HDM) of type II. Contributions contain top- and bottom-quark triangles, boxes 
and pentagon diagrams, i.e. the full mass dependence of the loop induced production. 
Interference effects between loops with bottom- and top-quarks as well as between CV- 
even and CP-odd couplings of the heavy quarks are fully taken into account. An option 
to use the large top mass approximation, which works well for intermediate Higgs boson 
masses (provided that the transverse momenta of the final state partons are smaller than 
the top-quark mass), is also implemented. 

Note that these gluon fusion processes processes require the use of the ggf lo executable 
rather than the vbfnlo one. 

If Higgs boson plus two jet production via gluon fusion is selected in vbfnlo.dat, 
i.e. PROCESS = 4100-4107., the following additional parameters can be adjusted in the 
ggflo.dat file: 



Q_LOOP : Input that sets how the Higgs boson coupling is determined. 

Effective theory in the large top-quark mass limit {rrit — ?■ oo). 

Coupling derived from top-quark loop. 

Coupling derived from bottom-quark loop. 

Coupling derived from both top- and bottom-quark loops. 



- Q_LOOP = 

- Q_LOOP = 1 

- Q_LOOP = 2 

- Q_LOOP = 3 

Default is set to 3. 



SUBPRQQ : Switch that determines whether the subprocesses with quark-quark initial 
state are included. Default is set to true. 

SUBPRQG : Switch that determines whether the subprocesses with quark-gluon initial 
state are included. Default is set to true. 

SUBPRGG : Switch that determines whether the subprocesses with gluon-gluon initial 
state are included. Default is set to true. 

HIGGS_MIX: Switch for Higgs mixing. If set to 0, there is no Higgs mixing, which is 
the default. If set to 1, the mixing is determined via the user-input variables: 

— CP_EVEN_MOD : Changes the strength of the CP-even coupling. Default is set 
to 1.0. 

— CP_ODD_MOD : Changes the strength of the CP-odd coupling. Default is set to 
1.0. 

using the Lagrangian 

-C Yukawa = Q (Ccvcn Vq + 1 C'odd 75 Vq) Q^ ■ (3) 
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with Cevcn = CP_EVEN_MOD and Codd = CP_ODD_MOD and 

Vb = yb = -rrib tan (3 

V 

Vt = yt = -mt cot (3 (4) 

V 

When working in the Standard Model, tan/3 is set to 1. Note that this option 
is not yet available when decays of the Higgs boson are included. The option 
HIGGS_MIX = 2 is only used when working in the MSSM, and incorporates mix- 
ing between all three neutral Higgs bosons according to the Z propagator matrix, 
as follows 

hi = Zih huSSM + ZiH -ffMSSM + ZiA ^MSSM (5) 

where i = HTYPE. 

4.4 susy.dat — parameters in the MSSM 

The file susy.dat is used to specify the supersymmetric parameters when working in 
the MSSM. This file is used either if MODEL = 2 in vbfnlo.dat, or if HTYPE = 1-3 in 
vbfnlo.dat. As described earlier, if MODEL = 2 the code will run in the MSSM, which 
means that the Higgs boson masses, widths and couplings will all be set to the MSSM 
values. If MODEL = 1 but HTYPE = 1-3, the code will run in the Standard Model, and 
only the Higgs boson mass will be affected. Consequently, this file affects processes 101- 
107 (Higgs boson plus 2 jet production via VBF), processes 110-117 (Higgs boson plus 
3 jet production via VBF) and processes 2100-2107 (production of a Higgs boson in 
association with a photon and 2 jets via VBF). In particular, the electroweak corrections 
to Higgs boson production via VBF are affected by the inputs of susy.dat. As stated 
in Section 13. 6[ the Higgs boson contributions to the production of three massive gauge 
bosons will be affected by susy.dat by fixing the Higgs boson mass and couplings, but 
the code will only provide a reasonable approximation to the full MSSM result in the 
decoupling region (i.e. when the Higgs is SM-like). 

• FEYNH.SWITCH : Determines whether FeynHiggs is used to calculate the MSSM 
Higgs boson sector. Default is false. When working in the MSSM, especially when 
including electroweak corrections, it is recommended that FeynHiggs is used, as 
the Higgs boson mixing and the renormalization scheme used in Vbfnlo are consis- 
tent with those used by FeynHiggs. If a SLHA file is used instead of FeynHiggs, 
inconsistencies may occur in the calculation of parameters. 

• SLHA_SWITCH: Determines whether the MSSM parameters are taken from a SLHA 
file. Default is true. If set to true the values in the SLHA input file are used 
instead of internal settings or calculationa^^l. 

• SLHA_FILE: Name of the SLHA input file. (An example file - spsla.slha - is 
provided, which is the default.) 

• BENCH : Various benchmark scenarios are pre-programmed in the code, as an alter- 
native to using a SLHA file. These are: 



^•"Note that a SLHA file can also be used as the input for FeynHiggs by setting both FEYNH_SWITCH 
and SLHA_SWITCH to true. 
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— : User-input scenario (see below for further inputs) 

— 1 : M^"^^ scenario jM] 

— 2 : No-mixing scenario J5T] 

— 3: Gluophobic scenario J5T] 

— 4: Small a^s scenario ^T\ 

— 5 : CPX scenario |52) 

— 10: SPSl£0 [M] 

— 11: SPSlb 

— 12-19: SPS2-9 

Default is 1 - the M™'^'' scenario. 

PROPLEVEL : Level at which Higgs propagator corrections are included: 



Effective Higgs-mixing angle used. 

Propagator factors included at leading order. 

Propagator factors included at leading order and loop level. 

Propagator factors included as an additional loop correction. 



Default is 1. These options are discussed in more detail on the Vbfnlo webpage, 
as well as in 1211. 



• DELMB_SWITCH : Switch determining whether or not to correct the bottom-quark 
Yukawa coupling. Default is false. 

If a SLHA file is not being used, and BENCH = 0-5, the following inputs are also needed. 

• TANB: Value of tan /3. Default value is 10. Note that for the SPS benchmarks (BENCH 
= 10-19), this value of tan /3 is not used. 

• MASSAO : Mass of CV-odd Higgs boson M^. This should be used when working in 
the MSSM with real parameters. Default value is 400 GeV. Note that for the SPS 
benchmarks (BENCH = 10-19), this value of Ma is not used. 

• MASSHP : Mass of charged Higgs boson Mh± ■ This should be used when working in 
the MSSM with complex parameters. Default value is -1 GeV. 

Note that it is standard to use as input the mass of the CV-odd Higgs boson, M^, when 
working in the MSSM with real parameters and the mass of the charged Higgs boson, 
Mh±, when working in the MSSM with complex parameter^. The mass that is not 
being used as input should be set to -1 GeV. 

If desired, the SUSY breaking parameters that define the scenario can be input by the 
user, by selecting BENCH = 0. Default values are those for SPSla. These parameters are 
(in the notation used by FeynHiggs): 



''The SPS points are defined with low-energy parameters, as given at 



http: //www, ippp.dur .ac .uk/~georg/sps/[ ). 

^"^This is because, when working with complex parameters, there is mixing between all three neutral 
Higgs bosons, and thus all neutral Higgs boson masses receive loop corrections. 
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• M3SQ etc.: The soft SUSY breaking parameters for the sfermion section. 

• AT etc.: Third generation trihnear couphng parameters. 

• MUE: Higgs boson mixing parameter. 

• M_l etc.: Gaugino mass parameters. 

4.5 Parameters for anomalous couplings 

Vbfnlo supports anomalous HVV couphngs, where V = W, Z, 7, in both the production 
and the decay of a Higgs boson in VBF-type reactions, i.e. for process IDs 100-107. The 
anomalous HVV couplings can be specified in the anom_HVV.dat input file. 

Anomalous triple and quartic gauge boson couplings are available for the VBF process 

pp — )■ W~^W~jj{j) |5l], the triple vector boson production processes WWZ, ZZW^ and 
W+W-W^ (process IDs 400-440) jSS] and for W-fj and WZj production I33ES]. The 
respective parameters are set in the input file anomV . dat[^. Note, however, that not all 
parameters in anomV . dat affect all processes. 

4.5.1 anom_HVV.dat — anomalous HVV couplings 

The file anom_HVV.dat controls the anomalous Higgs boson coupling parameters. It is 
used if the input ANOM_CPL in vbfnlo . dat is set to true. Among the anomalous coupling 
input parameters, the user can choose between three different parameterizations. 

The anomalous HVV couplings can be described by the following effective Lagrangian 
involving the dimension five operators 

HZZ HZZ HWW HWW 

C = ^^HZ^,Z^'' + ^^HZ^,Z^'' + ^^ HW+W!^'' + ^^ HW+W'^'' + 

/A5 2A5 A5 A5 

HZ^ HZ-y _ H77 H77 

^-^HZ^.A^" + ^if Z^,A^- + ^-^HA^,A^^ + ^-^HA^.A^" 

(6) 

where the subscript e or refers to the CV-even or CP-odd nature of the individual 
operators |56| . 

An alternative approach is to write the effective Lagrangian in terms of the operators 
Owwi Obbi Ow and Ob and their corresponding CP-odd operators according to Refs. [5] 



fw ^^ , Ib .^ , fww^^ , Ibb 
^^B + -jj 



C,E = 'ijOw + ^Ob + '-^Oww + ^Obb + CP-odd part (7) 

A^ Ki M Ki 



^^In previous versions of Vbfnlo, anomalous gauge boson couplings for WW production via VBF 
were input via the file ajiom_WW.dat. To ensure backwards compatibility, these old files can still be used 
(for this process only) by re-naming anom_WW.dat to anomV.dat. 
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The explicit form of the operators ii 



Ow = 


(D^0t)^/^-(/)^0) 


Ob = 


{D^<P^)B^''{D,<P) 


Oww = 


0^1^MvW^^"0 


Obb = 


(t>'^B^,B^''ct> 


Owww = 


Tr 


W^.W'PWp'' 


w,. 


= ^9T''W;, 


B^y 


= 


ig'YB^,, 



with 



(9) 

where g and g' are the SU(2) and U(l) gauge couphngs, T" are the SU(2) generators and 
Wfj_ij and B^ij are the field strength tensors. The CP-odd part of the Lagrangian has the 
same form, although only three parameters (denoted with a tilde) are needed. 

The different parameterizations, and the relationships between them, are discussed in 
more detail on the Vbfnlo webpage, where the explicit forms of the HVV couplings are 
given. 

1. A parameterization in terms of couplings in the effective Lagrangian approach given 
by Eq. O 

• PARAMETRl : Parameter that switches on the effective Lagrangian parameteri- 
zation of Eq. (E]). The default value is false. 

• LAMBDAS: Mass scale A5 in units of GeV with 480 GeV chosen as default. 

• G5E_HWW , G5E_HZZ , G5E_HGG , G5E_HGZ : Parameters that determine the cou- 
plings g^^^ of the CP-even dimension five operators. Their default values are 
set to 0. 

• G50_HWW, G50_HZZ, G50_HGG, G50_HGZ : Parameters that determine the cou- 
plings g^^^ of the CP-odd dimension five operators. Their default values are 
set to 0. 

2. The parameterization of the anomalous couplings by the L3 Collaboration as given 
in Ref. |59) . The parameters are d, ds-, Agf and Ak^, which are related to the 
coefficients /j/A^ of the effective Lagrangian of Eq. [7] in the following manner |60] : 

2 2 

"— A 2 JWW y "— .2 JWW^ 

_ m^ sin^6'^„ ~ _ rn^r sin^ 9^ 



'w 



Ak^ = K^ - 1 = ^ (/b + /w) , K^ = ^ /^ , 

Ao^ - 0^ - 1 - ^ ^— 
^9i -9i i - ^2 2 ■ 



(10) 
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Note that Owww is used only in the anomalous gauge boson couplings. 
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• PARAMETR2 : Parameter that switches on the above mentioned parameterization 
of Eqs. [ini The defauh value is false. 

• D_EVEN, DB.EVEN, DKGAM.EVEN, DG1Z_EVEN: Parameters that are the CV- 
even couphngs in this parameterization with defauh 0. 

• D_ODD, DB_ODD, KGAM.ODD: Parameters that are the CP-odd couphngs in this 
parameterization with defauh values equal to 0. 

• HVVl : Parameter that determines which anomalous HVV couplings are used 
for the run. For HVVl = 0, only the HZ'-y coupling, for HVVl = 1, only the if 77 
coupling, for HVVl = 2, only the HZZ coupling and for HVVl = 3, only the 
HWW coupling is used. If HVVl is set to 4, all possible anomalous couplings 
are used. This is also the default value. 

3. The parameterization of the anomalous couplings in terms of coefficients fi/h? of 
the effective Lagrangian in Eq. [71 

• PARAMETR3: Parameter that switches on the parameterization stated above. 
The default value is false. 

• FWW_EVEN, FBB.EVEN, FW.EVEN, FB_EVEN : Parameters that represent the co- 
efficients of the CP-even operators - i.e. /j/A^ - with default values equal to 
GeV-2. 

• FWW_ODD, FBB.ODD, FB.ODD: Parameters that are the coefficients of the CP- 
odd operators - i.e. /j/A^ - with default values GeV~^. 

• HVV2: Parameter that allows the user to choose which anomalous HVV cou- 
plings are used. For HVV2 = 0, only the HZ'^ coupling, for HVV2 = 1, only the 
if 77 coupling, for HVV2 = 2, only the HZZ coupling and for HVV2 = 3, only 
the HWW coupling is used. If set to 4 all possible anomalous couplings are 
used. The default value is 4. 

Moreover, for all parameterizations, two different form factors can be chosen as described 
in Refs. |^l^. They model effective, momentum dependent HVV vertices, motivated 
by new physics entering with a large scale A at the loop level. 

A^ A2 

F2 = -2K'Co{qlql{qi + q2fA'). (12) 

Here the qi are the momenta of the vector bosons and Cq is the scalar one-loop three 
point function in the notation of Ref. |61] . 



FORMFACTOR : Parameter that switches on the above form factor. The default value 
is set to false. 

MASS_SCALE: Characteristic mass scale. A, of new physics in units of GeV. The 
defauh value is 200 GeV. 

FFAC : Parameter that is used to select one particular form factor out of Eqs. (ITT]) 
and f lT2|) . If FFAC = 1, the form factor Fi is used for the parameterization. FFAC = 
2 selects i^2, which is also the default value. 
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Finally, the following two parameters can be used to rescale the SM HVV couplings. 

• TREEFAC : Parameter that multiplies the HVV tensor present in the SM Lagrangian. 
Default is 1. 

• LOOPFAC: Parameter that multiplies the HZ'-f and H-f-f vertices induced by SM 
loops. The default is chosen to be 1. 

4.5.2 anomV.dat — anomalous triple and quartic gauge boson couplings 

The anomalous triple and quartic gauge boson couplings can be set in the file anomV.dat. 
They are parameterized using an effective Lagrangian, as described in Refs. [ ^ 1 ^ 1^ 

^cfT = ^or\ (13) 

where n + 4 signifies the dimension of the operator Oi. Vbfnlo defines the anomalous 
gauge couplings in terms of the coefficients /j/A" of the dimension 6 and (for triboson 
production) dimension 8 operators (a full list of implemented operators for triboson pro- 
duction can be found in Ref. |55|). 

A common alternative parameterization (which Vbfnlo can also use as input) of the 
trilinear couplings WW'-f and WWZ uses the following effective Lagrangians: 



^WWy — — ^6 

for the anomalous WW'-f vertex, and 
^wwz = — ie cot 9. 



wl^WA" - wlA^w"'' + K^wlw^Ff"" + ^w^^^wi;f' 



^w 



(14) 



9i {Wl^Z' - WlZ.W^"") + KzWlW.Z^'' + ^Wl^Wl^Z' 



m^w 



(15) 



for the anomalous WW Z vertex. It is customary to rephrase the electroweak modifications 
around the SM Lagrangian in terms of new quantities, 

(A(7f , Aavz, Afi:^) = (^f , /tz, ft:^) " 1- (16) 

These quantities are related to the coefficients of the dimension 6 operators O^y, Ob 
and Ovi/vyvK in the following fashion. 



^9i 



z _ JW^^^Z 



fwMl 
2 

1\kz = COS ow^rw ~ sm "wt^ 



Ak^ 



A2 2 ' 

fw . 2n fB\Ml 

— - sm Owj^ 
fw , Ibx M^ 



A2 A2/ 2 



\ _ \ _ 3Mj^e fwww /.„x 

^' - ^^-2sin2^H. A^ ■ ^ ^ 

In order to include anomalous vector boson couplings, the parameter ANOM_CPL must 
be switched to true in vbfnlo.dat. The parameters described above are input via the 
file anomV.dat: 
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TRIANOM: Switches between parameterizations of the anomalous WW'y and WWZ 
couphngs. TRIANOM = 1 uses the coefficients of the dimension 6 operator^j as input: 

— FWWW: Coefficient of the operator Owww, i-e. fwww/-^"^- Defauh is set to 
9.19 X 10-6 GeV-2. 

— FW: Coefficient of the operator Ow, i-e- fw I ^ ■ Defauh is set to —1.44 x 
10-5 GeV-2. 

— FB: Coefficient of the operator Ob^ i-e. JbI^- Defauh is set to 3.83 x 
10-5 GeV-2. 

TRIANOM = 2 uses the ahernative parameterization of Eqs. \\M and [TSl as input: 

- LAMBDAO: The quantity \{= \z). Defauh is set to 0.038. 

- ZDELTAKAPPAO : The quantity Ak^- Defauh is set to -0.082. 

- ZDELTAGO: The quantity A^ff. Defauh is set to -0.060. 

- ADELTAKAPPAO : The quantity Ak^. Defauh is set to 0.077. 

— Note that, as can be seen from Eqs. [T71 the quantities Ak^, Ak^ and A(7f are 



not independent, but obey the relation 



sin^ d\ 



. A 7 sin (7W . /, ^N 

Akz = A(?f ^Aft:^. 18 

COS^ t^VF 

If one of these quantities is zero, it will be set by Vbfnlo to be consistent 
with the other values. If the input values are inconsistent, Ak^ will be reset to 
give the correct relation. 

Default is TRIANOM = 1, and the default values are approximately the central values 
in the three-parameter fit given in Ref. [6i 



• FEW, FDW, FWW, FEB: The coefficients of the remaining CP-even dimension-6 op- 
erators, i.e. /j/A^. Note that these are only implemented for W pair production via 
vector boson fusion. Default value is GeV"^. 

• FWWWt, FWt, FEt, FEWt, FDWt , FWWt , FEBt : The coefficients of the CP-odd dimension- 
6 operators, i.e. /j/A^. Note that these are only implemented for W pair production 

via vector boson fusion. Default value is GeV"^. 

• FSO, FSl, FMO - FM7, FTO - FT2, FT5 - FT7 : Parameters that give the values 
of the coefficients of the dimension-8 operators, i.e. /i/A^. The default values for 
these parameters are GeV"^. Note that these are only relevant for triboson pro- 
duction. 

In addition, a form factor can be applied: 

i^=(l + X^) , (19) 

where A is the scale of new physics. 

^^The explicit form of the operators is given in Eq. [51 
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• FORMFAC : Switch determining whether the above form factor F is included. Default 
is set to false. 

• FFMASSSCALE: Mass scale A. Default is set to 2000 GeV. 

• FFEXP : The exponent p. Default is set to 2. 

Individual form factor mass scales and exponents can be set for the trilinear WWZ and 
WW'j couplings, using either of the parametrizations above. If chosen, these values 
overwrite the "universal" mass scale and exponent (FFMASSSCALE and FFEXP) set above 
for the selected parameters. 

• FORMFAC_IND : Switch determining whether individual or universal form factors are 
used. Default is set to false - universal form factors are used. 

If TRIANOM = 1 then 

• MASS_SCALE_FWWW: Mass scale A for coefficient fwww- Default is 2000 GeV. 

• FFEXP_FWWW: Exponent p for coefficient fwww- Default is 2. 

• MASS_SCALE_FW: Mass scale A for coefficient fw Defauh is 2000 GeV. 

• FFEXP_FW: Exponent p for coefficient fw Default is 2. 

• MASS_SCALE_FB: Mass scale A for coefficient fs- Defauh is 2000 GeV. 

• FFEXP_FB: Exponent p for coefficient fs- Default is 2. 
If TRIANOM = 2 then 

• MASS_SCALE_KAPPA : Mass scale A for parameter An. Default is 2000 GeV. 

• FFEXP_KAPPA: Exponent p for parameter An. Default is 2. 

• MASS_SCALE_LAMBDA: Mass scale A for parameter A. Default is 2000 GeV. 

• FFEXP_LAMBDA : Exponent p for parameter A. Default is 2. 

• MASS_SCALE_G: Mass scale A for parameter Ag. Default is 2000 GeV. 

• FFEXP_G: Exponent p for parameter Ag. Default is 2. 



Finally, an overall factor scheme |S1] can be employed for W^W production via VBF, 
using the switch: 

• OVS: Switch for overall factor scheme. Default is false. 
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4.6 kk_input.dat — parameters for Higgsless models 

Vbfnlo allows the calculation of the weak boson fusion processes WWjj, WZjj and 
ZZjj in the Warped Higgsless scenario p5|IB5] at LO and NLO QCD level (see, for ex- 
ample, Ref. |67] for a phenomenological application). The triple vector-boson processes 
WWW, WW Z and W ZZ can be calculated in both this model and a Three-Site Higgs- 
less Model |68| . The model parameters can be generated by Vbfnlo via the input file 
kk_input.dat for a choice of the relevant five dimensional gauge parameters. The input 
parameters are: 

• VBFCALC.SWITCH: Switch that determines whether Vbfnlo should calculate the 
parameters needed by the model. Default is true. 

• MHLM_SWITCH: Switch that determines whether the Three-Site Higgsless Model [68] 
(true) or the Warped Higgsless Model |65] (false) should be used. Default is false 
(Warped Higgsless Model). 

• The following parameters are used only when working in the Warped Higgsless 
Model (MHLM_SWITCH = false): 

— RDWN : Location of the UV brane for the generation of the model parameters in 
the Warped Higgsless Model. Default is i? = 9.75 x 10~^, which amounts to the 
Kaluza-Klein excitations having masses of m^y^ = 700 GeV, mz^ = 695 GeV, 
and rriz' = 718 GeV. Smaller values of R result in a heavier Kaluza-Klein 
spectrum. 

— KKMAXW: The maximum number of Kaluza-Klein W^ states to be included on 
top of the Standard Model W"^ bosons, which correspond to W^^^q- ^^^ states 
A; > 3 are phenomenologically irrelevant. Default is 1. 

— KKMAXZ : The maximum number of Kaluza-Klein Z^ states to be included on 
top of the Standard Model Z boson, which corresponds to Zk=o. All states 
k > 3 are phenomenologically irrelevant. Default is 1. 

— KKMAXG : The maximum number of Kaluza-Klein Z(, bosons that are excitations 
of the Standard Model photon Z'^^^q. States k > 2 are phenomenologically 
irrelevant. Default is 1. 

• WPMASS: W^fc^i mass in the Three-Site Higgsless Model (MHLM.SWITCH = true). De- 
fault is 500 GeV. 

The explicit breaking of higher dimensional gauge invariance is balanced according to 
the description of Ref. 1^1], where more details regarding the model and its implementation 
can be found. Information on the widths and the sum rules relating the various gauge 
boson couplings j70] are written to the file kkcheck.dat. 

Vbfnlo generates the text file kk_coupl_inp.dat, which documents the calculated 
model parameters, i.e. Kaluza-Klein gauge boson masses and couplings of the specified 
input parameters. This file can also be used as an input file for advanced users who 
want to run the code with their own set of parameters - e.g. for simulation of more 
general technicolor-type scenarios. To that end, select VBFCALC_SW1TCH = false in the 
file kk_input.dat. Vbfnlo will then calculate the gauge boson widths on the basis of 
these parameters by the decay to the lower lying states. The inputs in kk_coupl_inp . dat 
are as follows: 
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KKMAXW: The niaxiniuin number of Kaluza-Klein Wf^ states included on top of the 
Standard Model W^ bosons. 

KKMAXZ : The maximum number of Kaluza-Klein Zk states included on top of the 
Standard Model Z boson. 

KKMAXG : The maximum number of Kaluza-Klein Z'j. bosons that are excitations of 
the Standard Model photon Z'i^^q. 

Kaluza-Klein masses: 

- KKMASSW: List of masses of the Kaluza-Klein W^ states. 
List of masses of the Kaluza-Klein Zu states. 



KKMASSZ 
KKMASSG 



List of masses of the Kaluza-Klein Z'^. states. 



Couplings 



CPL_Wx-Wy-ZX : List of couplings between Kaluza-Klein states: W^ — Wy — Zq 

to Wt - W^ - Zkkmxz. 

CPL_Wx-W|/-GX : List of couplings between Kaluza-Klein states: W^ — W^ — Z'q 
to W^ - W^ - Z' 

CPL_Wa;-Wy-W2;-WX : List of couplings between Kaluza-Klein states: W^ — W^ — 

\Y± _ \Y± to W^ - W^ - W^ - W^ 

''^ Z ''^ ^^ ''^ X ''^ y "^ Z '" KKMAXW • 

CPL_Wx-Wi/-Z2;-ZX : List of couplings between Kaluza-Klein states: W^ — W^ — 
Z, - Zo to Wt - W^ -Z,- Z^xz. 

CPL_Wx-Wi/-Gz-GX : List of couplings between Kaluza-Klein states: W^ — W^ — 
Z'^ - Z', to wt - Wf -Z'^- Z4_. 

CPL_Wx-Wi/-G2;-ZX : List of couplings between Kaluza-Klein states: W^ — W^ — 
Zi - Zn to W^ - Wf -Z[- Z. 



JQ \.U yV — yVy — ^z~ -^KKMAXZ- 



The default values in kk_coupl_inp.dat are those produced when Vbfnlo is used to 
calculate the couplings and masses using the default values of kk_ input . dat as described 
above. 

4.7 histograms.dat — parameters for histogram ranges 

The file histograms.dat allows the user to set the range of the x-axis of the pro- 
duced histograms. The range is input in the format: xmin xmax. The following de- 
scribes those histograms which are already implemented in Vbfnlo. By altering the file 
utilities/histograms .F, however, it should be easy for the user to add new histograms 
- Vbfnlo will automatically read-in the range of each created histogram. 

• HIST_ID1 : Range for pt of tagging jets. Default range is to 250 GeV. 

• HIST_ID2: Range for p^^ of tagging jet with higher p'p. Default range is to 500 GeV. 

• HIST_ID3 : Range for pj- of tagging jet with lower p^- Default range is to 200 GeV. 

• HIST_ID4: Range for pseudorapidity of tagging jets. Default range is -5 to 5. 
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• HIST_ID5: Range for pseudorapidity of tagging jet with higher pt- Default range 
is -5 to 5. 

• HIST_ID6 : Range for pseudorapidity of tagging jet with lower pt- Default range is 
-5 to 5. 

• HIST_ID7 : Range for p'f^'' of leptons. Default range is to 500 GeV. 

• HIST_ID8: Range for p^?'"^ of leptons. Default range is to 500 GeV. 

• HIST_ID9 : Range for ?7max of leptons. Default range is to 5. 

• HIST_ID10: Range for ?7min of leptons. Default range is to 5. 

• HIST_ID11 : Range for azimuthal angle of tagging jets. Default range is -180 to 180. 

Vbfnlo can also produce 2D histograms. In this case, both the x and y ranges can 
be set using histograms. dat in the format xmin xmax ymin ymax. 

• 2DHIST_ID1 : Range for 2D histogram of (PS/{drijjd'mjj). Default ranges are to 6 
for the f] (x) axis and to 800 GeV for the rrijj [y) axis. 
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5 Changes 

The current version of Vbfnlo - 2.5.0 - has been altered in such a way that some results 
will differ from previous versions. These changes are described briefly here, and in more de- 



tail on the Vbfnlo website, http://www-itp.particle.uni-karlsruhe.de/vbfnlo/ 



5.1 EWSCHEME 

The implementation of options EWSCHEME = 1 and 4, which are described in Section 14.1.2^ 
has been altered. The new implementation is hopefully more transparent, and is described 
in this manual. 

5.2 Gluon fusion 

In the current version of Vbfnlo, the gluon fusion cross section has also been affected by 
changes to the code. In the new version, the bottom-quark mass Mi,{Mh) is used through- 
out the calculation. New, more stable, tensor routines for the boxes and pentagons have 
been implemented, reducing the number of omitted points with bad Gram Determinants. 
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6 Checks 

Extensive checks for the LO and the real emission aniphtudes as well as for the total 
LO cross sections have been performed for all processes implemented in Vbfnlo. Born 
amplitudes and real emission diagrams have been compared with the fully automatically 
generated results provided by MadGraph jTT]. Complete agreement has been found 
in each case. Moreover, total LO cross sections with a minimal set of cuts agree with 
the respective results obtained by MadEveniH [T21[73], Sherp/@ |I1] and/or Helac- 
PHEGAq^ [75H77] . a completely automatic parton level event generator based on Dyson- 
Schwinger recursive equations. 

LHA event files for the LO processes have been tested with HerwigH — hcj [ZH], a 
general purpose Monte Carlo event generator for the simulation of hard lepton-lepton 
and hadron-hadron collisions. 

As a final and very important test, comparisons with already published results have 
been made. The NLO results for Higgs boson production via VBF agree with those 
produced by the code Hawko. In Ref. |79| . a tuned comparison of LO and NLO QCD 
results for Higgs boson production via vector boson fusion at the LHC has been performed. 
Three different calculations have been cross checked: Vbfnlo, the results of Refs. JMIEI], 
and the VV2H prograno. For the dominant t- and w-channel contributions which are 
implemented in Vbfnlo, good agreement has been found. For the triboson processes a 
comparison for the production of on-shell gauge bosons without leptonic decays has been 
performed with the results presented in Ref. [82]. Additionally, the tree-level results for 
W^'-ff have been compared to the results with an on-shell W boson of Ref. [83] • Again, 
in both cases good agreement has been found. Triple photon production has been tested 
against FeynArts, FormCalc and HadCalc [3llll|8l]. Results for the CP-odd and 
CP-even Higgs boson production via gluon fusion have been tested against FeynArts 
and FormCalc I85ll86l. 
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It^^T* 



http : //madgraph . hep . uiuc . edu/ 
http : //www . sherpa-mc . de/J 



http: //helac-phegas .web. cam. ch/he lac-phegas/| 
http : //pro j ects . hepf orge . org/herwig/ 
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http://oinnibus.uni-freiburg.de/~sd565/prograins/hawk/hawk.html 



http: //people .web.psi . ch/spira/vv2h/ 
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7 Outlook 

Additional processes will become available in the near future and will be included in the 
code version on the Vbfnlo website. 
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